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Abstract 

We suggest an unconventional view of the origin of most cosmic rays (CRs) in the 
universe. We propose that nearly every accelerated CR was part of the parallel current 
that maintains all force-free (f-f) magnetic fields. Charged particles are accelerated by 
the Eu (to the magnetic filed B) produced by reconnection. The inferred total energy 
in extra-galactic cosmic rays is ~ 10 60 ergs per galaxy spacing volume, provided 
that acceleration mechanisms assumed do not preferentially only accelerate ultra high 
energy cosmic rays (UHECRs). This total energy is quite large, about x 10 5 the parent 
galactic CR or magnetic energy. We argue that the formation energy of supermassive 
black holes (SMBHs) at galaxy centers, ~ 10 62 ergs, becomes the only feasible source. 

We propose an efficient dynamo process which converts gravitational free energy 
into magnetic energy in an accretion disk around a SMBH. Aided by Keplerian wind- 
ing, this dynamo converts a poloidal seed field into f-f fields, which are transported into 
the general inter-galactic medium (IGM) eventually. This magnetic energy must also 
have been efficiently converted into particle energies, as evidenced by the radiation 
from energetic particles. In this view CRs of the IGM are the result of the continuing 
dissipation, in a Hubble time, of this free energy, by acceleration in situ within the f-f 
fields confined within the super-galactic walls and filaments of large scale structures. 
In addition, UHECRs are diffusively lost to the galactic voids before the GZK attenua- 
tion time, ~ 10 8 years. Similarly within the galaxy we expect that the winding by the 
disk rotation of the galaxy, by the rotation energy of magnetized neutron stars, and 
by the Keplerian winding of star formation disks are efficient sources of f-f magnetic 
field energy and hence the sources of galactic CR acceleration. 



1 Introduction 



In this view of the origin of CRs, the source of the necessary free energy and the form 
that it takes, force-free magnetic fields, is the organizational principle in determining the 
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mechanism of CR acceleration. Here not only the CR acceleration mechanism is different 
from traditional views, but also the implied and necessary strength and origin of the 
magnetic fields in the general IGM. There are literally hundreds of topics and issues to be 
addressed, both on the positive and the negative side of this view. We will address a few 
of these here, but recognize that attempting to reverse the conclusions of so many people 
over so many years is well beyond this single article. However, we emphasize again that 
the available free energy for CR acceleration is the organizing principle of this view of CR 
origin. 

There are three primary issues: 

1) The likely total energy of extra galactic CRs (including radio lobes) compared to the 
likely sources of this energy. 

2) The astrophysical circumstances for the generation of this energy as f-f magnetic field 
energy within the galaxy and within the meta galaxy. 

3) The scaling in time and space of E\\ reconnection acceleration leading to the CR spec- 
trum and energy upper limit. 



2 Total Energy in Extra Galactic and Galactic CRs 
2.1 CR Spectra and Implied Energies 

There are two circumstances where the evidence for extra galactic cosmic rays appears to 
be beyond doubt, UHECRs and giant radio lobes. The first is because the particle energy 
is great enough, beyond the "ankle" or E > 10 18 eV, that confinement by the magnetic 
fields of the Galaxy is unlikely. Here the spectral index of UHECRs returns to nearly the 
same slope as at lower energies, below the "knee" at E < 10 15 eV. A likely explanation 
of the change in slopes, T = —2.7, —3.0, and — 2.6, for the energy regions E < 10 15 , 
10 15 < E < 10 18 , and > 10 18 eV, is that at lower energies, proton escape from the galaxy 
at E < 10 15 eV, progressively higher A/ < Z > nuclei escape < 10 15 < E < 10 18 eV 
up to iron at E ~ 10 18 eV and extra galactic protons above this energy. The expected 
enrichment at < 10 15 < E < 10 18 eV has been observed (Abu-Zayyad et al. 2001). 
This spectrum is shown in Figure 1 with various extrapolated slopes superimposed from 
E < 10 18 down to rripC 1 . Normalizing at 10 19 eV and extrapolating the spectrum back 
to E = nipC 2 using a slope of T = —2.0 gives a minimal estimate of the total energy 
in extra galactic CRs. The result is that the local specific energy content in the extra 
galactic component in the meta galaxy, e mg , is e mg ~ 2 x 10~ 19 erg/cm 3 , as compared to 
€ ga i ~ 10~ 12 erg/cm 3 for the energy content of the galactic cosmic ray spectrum. Since the 
ratio of the two volumes is ~ 4 Mpc 3 / ~ 300 kpc 3 ~ 10 7 , both total CR energies, inside 
and outside the galaxy, would be about equal. However, such a flat spectrum as V = —2.0 
is quite optimistic for any stochastic acceleration process. It implies that the total energy 
in accelerated particles per logarithmic interval in energy is a constant so that the integral 
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or total energy diverges with the upper energy limit. The back reaction from such an 
efficient acceleration mechanism would necessarily alter the source of the energy leading 
to a self limiting spectrum. Values of T = describe a modern research accelerator where 
losses during acceleration are near zero and the spectrum is a delta function in energy. The 
astrophysically reasonable spectrum is for T < —2.0. Acceleration to such high energies 
has confounded astrophysicists for nearly a century. The slope, T = —2.7, is constant over 
6 decades, for CRs of energy Eqr < 10 15 eV within the galaxy and near constant extra 
galactic over 2 decades with slope T ~ —2.6 for Eqr > 3 x 10 17 eV. The constant slopes 
of closely the same value strongly imply a single scale independent mechanism, so that for 
extra galactic CRs we expect the slope F ~ —2.6 to extend as a constant to low energy 
the same as in the galaxy. The small difference in slopes is likely explained by the small 
(over 6 decades) additional nuclear scattering spallation loss, ~ 50% in 10 7 years within 
the galaxy. 

With this assumption in Fig. 1 we have drawn extrapolated spectra backwards at various 
slopes from the measured values at high energy. First we note that at the high energy 
where the galactic magnetic field confinement should cease to modify the spectrum, the 
"ankle" or E > 10 18 eV, defining CRs in the meta galaxy, the ratio of energy densities 
becomes e mg /e ga i ~ 1/30. This value is comfortably close to what we would expect if 
galactic and extra galactic CRs reach equipartition in energy with what we believe are 
the magnetic energies within and external to the galaxy. We will discuss the origin and 
evidence for the extra galactic fields later in more detail. 



2.2 Radio Lobes 



A second indicator of extra galactic cosmic rays is extra galactic radio lobes. Giant 
radio lobes are a strong signature of extra galactic in-situ particle acceleration. Since 
radio lobes are detected by their radio emission and since strong polarized correlated 
emission is observed over distances of up to Mpc, only synchrotron emission by relativistic 
electrons in a magnetic field makes a sensible explanation. However, it has been known 
since Burbidge, (1956) that a minimum in the total energy required to create the radio 
lobe emission, (1) wave length, (2) luminosity, and (3) dimension, occurs for specific 
values of magnetic field, electron number and electron energy. These total energies as 
reported for up to 70 such giant radio lobes in Kronberg et al. (2001), are immense, up 
to 10% of the MBH rest mass energy. In that analysis the minimum energy is calculated 
as if the electrons are accompanied during acceleration by xlOO the same number of 
protons at a given energy, just as observed in the Galaxy. However, to be conservative, 
not knowing the acceleration mechanism in detail, we assumed, as did others, that the 
spectrum commences at the minimum energy necessary of the emitting electrons, 7 = 
E e /m e c 2 ~ 3x 10 4 , or E p ~ 3x 10 13 eV with negligible low energy particles, E p / (m p c 2 ) ~ 1. 
If we believe that it is more likely that only one acceleration mechanism prevails for all 
energies because of the constancy of T, or equivalently is most efficient, then extrapolating 
down to the total energy maximum of the spectrum, m p c 2 corresponds to an increase of 
the proton total energy of xl0( 4 - 5 ) r 2 = 10 2 - 47 = x300. However, the total minimum 
energy of the lobe increases only as the 4/7's power of the energy of any one component so 
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Figure 1: The many measurements of cosmic rays at high energy as compiled by Gaisser and 
Stanev (2000). The data are plotted as E 2J dN/dE vs E and thus with the slope V = -2.7 
removed. The change in slope to T ~ —3.0 in the energy interval E < 10 15 < E < 10 18 > 10 18 
eV is evident between the " knee" and " ankle" . This change in slope is interpreted as a transition 
from galactic confinement to extra galactic CRs. The width of the transition region is interpreted 
as progressive loss of nuclei of progressively higher < A > /Z with the highest being iron nuclei. 
Particles above 10 18 eV are interpreted as extra galactic UHECRs, primarily protons. Several lines 
are drawn at different slopes tangent to the UHECR spectrum and extrapolated back in energy to 
Ecr = rripC 2 where the maximum in total CR energy resides. The slope T = —2.7 is the galactic 
slope and results in an energy density of extra galactic CRs is ~ 7% of the galactic value. The 
slope for CRs in the IGM for the same acceleration mechanism should be slightly flatter because 
of the lack of spallation. Spallation attenuates galactic CRs by ~ 1/2 in 6 decades in energy and 
so we estimate that without spallation T ~ —2.6. this results in the energy density of ~ 0.005 of 
galactic CRs. 

that minimizing, but including the low energy extrapolation, increases the minimum total 
energy by the factor x 10^ 2 ' 47 ^ <4/7) = 10 1 ' 4 = x25 when the new minimum is sought. If 
one chooses not to accelerate the protons, the minimum energy is reduced by ~ 10. Thus 
depending upon assumptions of non-standard acceleration, the giant radio lobes contain 
CRs up to a hundred times the flux of CRs in the IGM or contain only slightly less energy 
than the maximum possible back hole formation free energy. We conclude that the extra 
galactic CRs require so much energy, that only the free energy of black hole formation 
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is a feasible source and that furthermore the radio lobes are a signature of this energy 
and acceleration mechanism. Finally this picture of immense high energy electron fluxes 
within radio lobes has recently been confirmed by Chandra observations of x-rays from 
high 7 ~ 10 3 electrons Compton scattering cosmic background photons (Harris 2003). We 
have concluded separately, Kronberg et al. 2004, that these radio lobes must be regions 
of immense f-f fields in quasi magneto hydrostatic equilibrium. 

2.3 Our Basic Model 

The magnetic field within the Galaxy is well recognized to be B ga \ ~ 5/xG, but we will claim 
in this paper the equally unconventional large magnetic field in the meta galaxy, B mga i ~ 
1/iG. We will discuss further this estimated large value of the magnetic field in the meta 
galaxy later, but here note that this ratio of magnetic energy densities is consistent with the 
expected infall pressures, presumably a confinement pressure, of the IGM accreting onto 
the galaxy and the infall pressure from the matter in the voids accreting onto the galaxy 
walls or filaments. At this accretion rate the masses of both the galaxy and filaments 
increase by ~ 50% in a Hubble time. The extra galactic field of B mga i ~ lfiG is also 
consistent with the inferred radio lobe magnetic fluxes derived from minimum energy 
derivations, rotation measure limits, and the theoretical magnetic flux production from 
the SMBH accretion disk dynamo. When we slightly increase the extra galactic index to 
r = —2.6 from the value T = —2.7 within the galaxy, the total energy density in extra 
galactic CRs is then reduced to e mg ai/e ga i ~ 1/300. This is sufficient such that the total 
energy in extra galactic CRs per galaxy spacing volume of ~ 4 Mpc 3 becomes (1/300) x 10 62 
ergs or (1/300) x 10 8 Mqc 2 , where 10 8 M & c 2 is the energy available in the accretion disks 
that form the SMBHs. Thus the local, (to the wall or filament) extra galactic CRs may 
be replenished ~ 100 times in a Hubble time. This replenishment time, ~ 10 8 years, is 
close to the estimated loss time by random walk using a semi-coherent intergalactic field 
of a galactic spacing (GS) distance, das — 2 Mpc and a filament thickness of 5 GS or 
5 galaxies thick or 10 Mpc. The random walk distance to escape the filament becomes 
^filament/ ^GS — 30 Mpc. This leads to a typical loss time of the CRs to the voids of 
ti oss ~ D"j ilament /cd.GS — 10 8 years. This loss or replenishment time ensures that roughly 
half of the UHECRs that are subject to the GZK loss will have been lost to the voids before 
detection at Earth provided the total energy available for continuing acceleration during 
a Hubble time is ~ 10 62 ergs. This loss time is also the maximum likely, upper limiting 
acceleration time since the f-f fields are primarily confined to the filaments. The likely 
upper energy limit of Eh acceleration for a particle that remains in a reconnection flux tube 
for the entire time will be E max ~ B[x x ti oss c = 3 x 10 22 eV where [i = 300, the conversion 
from gauss to volts at a velocity c. This limit excludes any possible reconnection current 
filamentation that may locally increase the acceleration field, En , but the Larmor radius 
within the filament and dependence upon its size and geometry in 10~ 6 /i G field may limit 
the maximum energy to an order of magnitude less, ~ 3 x 10 21 eV. In this picture there 
are no local sources, just a space-filling acceleration and a possible bias in flux opposite 
to the direction of the nearest void. 

In this picture CRs are accelerated by the same mechanism inside the Galaxy as well 
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as extra galactic by a space-filling reconnection of f-f fields; they are closely isotropic, 
mostly not attenuated by cosmic background photons, diffuse in an equipartition field 
both inside and outside the galaxy, where the flux and field are consistent with long time 
mass accretion rate estimates and are consistent with an energy source of SMBH accretion 
disks, extra galactic and similar sources within the Galaxy. We therefore discuss next the 
basis of f-f fields in the astrophysical environment and the more limited possibilities of 
magnetic field energy generation from the free energy of SMBH formation. 



3 Force-Free Fields in Astrophysics 



The anzatz of this paper is that once an accretion disk has been "seeded" with magnetic 
flux of sufficient magnitude, then subsequent accretion will wind up this flux, polloidal, into 
a force-free field, a combination of both toroidal and polloidal, whose total energy is that 
of the gravitational energy released in accretion. The seed flux is produced by the dynamo 
action occurring either in the MBH accretion disk, or in stars (Colgate, Li & Pariev 2001). 
The physics of the a u dynamo is beyond the subject of this paper, but without such a 
dynamo, it is most unlikely that primordial effects such as the Bierman battery or early 
phase transitions can supply the seed field necessary to transfer the angular momentum 
of accretion to tension in magnetic field. The energy contribution of the dynamo is small 
just as the exciter field of any commercial electric generator is a trivial fraction of the 
power produced, but without an exciter field no power is produced. The subsequent flow 
of power by the winding of the polloidal field within a conducting medium, external to 
an accretion disk, a Poynting flux, produces the helical f-f field that we associate with 
the "jet" of MBHs or jet identified with star formation (Lovelace , Li et al. 2001). Fig.2 
shows a composite of how such an accretion disk, dynamo, and f-f helix are formed. The 
f-f magnetic helix is a minimum energy configuration where the major magnetic stress is 
maintained by the tension in the field itself. The f-f helix is then self illuminating by the 
electrons accelerated by reconnection and appears as a jet in shape only. The relativistic 
velocities are then a combination of the relativistic synchrotron emitting electrons and the 
phase velocity of the reconnection instabilities. 

The number density of the current carriers, n e necessary to carry the current, Jm, is 
trivially small compared to the surrounding IGM, and so the mass carried by such a jet 
is also trivially small. The "beta" of relativistic current carriers Vdrift — c at radii ~ 3r g , 
r g = MG/c 2 , the gravitational radius, is (3 = n e m e c 2 / (B 2 / 8tt) ~ 10" 15 . Such a f-f field 
jet is then a self illuminated (in x-rays, optical, and radio) f-f field configuration. It is 
illuminated in synchrotron radiation by electrons accelerated the same as in the radio lobes. 
A similar jet is formed in star formation or by the winding of flux by either neutron stars, 
as we believe in the Crab nebula or by the galaxy itself (not yet observed but predicted). 
In any case the jet subsequently morphologically transforms by tearing mode reconnection 
into the radio lobes and finally into a galactic and inter galactic space-filling f-f magnetic 
flux and therefore a source of free energy. It is the continuing reconnection of this space- 
filling f-f magnetic flux that we associate with the production of a near equipartition flux 
of galactic and extra galactic CRs. Because of the general lack of familiarity with f-f fields 
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in astrophysics, we give a laboratory example of the physical steps necessary to produce 
such fields and point out the analogy with accretion disks. We conclude this section with 
a description of the evolution of this flux within cosmic structure calculations. 

3.1 A Laboratory Example of the Generation of Force Free Fields 

The magnetic field inside a current carrying coil is orthogonal to the boundary current and 
the force of the field on the coil winding is J x B = B 2 /8ir, which we will call a "force- 
bounded" magnetic field. In particular J • B ~ everywhere except possibly the coil 
current leads. If now two conducting metal disks are located at either end of the coil and 
are threaded by the predominantly axial field, the field will penetrate the conducting metal 
disks in a finite time and the field will return to its original primarily axial configuration. 
If now the disks are rotated relative to each other, an electric field will develop both 
radial and axial, but nothing further will happen because the medium between them, air, 
is an insulator. If now the air is replaced by a conducting fluid medium, e.g., plasma or 
liquid metal and provided the conductivity, disk rotation velocity, and dimension are great 
enough, i.e., the magnetic Reynolds number, Rm = vL/rj >> 1, a current, I z will flow 
both plus and minus axially at inner and outer radius. The axial currents are connected 
to radial currents, Ir within the disks. The axial current will produce an azimuthal field, 
bounded by the plus and minus axial currents as well as by the plus and minus radial 
currents within the disk. The vector sum of the axial and azimuthal fields B z + B^ 
changes the topology of the field from initially axial to helical. Furthermore the radial 
current in the disks produce a torque, R ■ Jr, x B z , which, times the winding of the disks, 
performs the work that increases the magnetic energy. In addition this field will exert a 
radial force on the conducting fluid depending upon its confinement, a "pinch". The same 
plasma configuration, a "stabilized pinch" is produced in the laboratory by replacing the 
rotation of the disks by a voltage (capacitors) between the disks and a low pressure plasma 
between them. The resulting current between the disks (electrodes) produces the current 
of the resulting helical, f-f fields. If the pressure of the conducting plasma is small or some 
of the incompressible liquid metal can escape either because of instabilities or through the 
disk wall, then the fluid can not exert a force on the field due to the induced azimuthal 
field or original axial field. The resulting configuration is called a force-free field where 
without force, J x B = and therefore J = AB. In general when the minimum energy 
state is reached of such a f-f field, subject to various boundary conditions of flux and 
helicity conservation, the field strength decreases with radius B ~oc 1/R so that a much 
weaker outer boundary field exists which must be supported at some conducting boundary 
with a weak J x B = V(B 2 /8-7r) = VP force. This boundary pressure in the astrophysical 
case is the pressure of the IGM or ISM. However, the main point is that the differential 
rotation of the disks has added free energy to the field by the work done by twisting the 
field, (equivalently supplied by the capacitor in laboratory plasma experiments) thereby 
increasing both components of the field, (B 2 + B^)/8it. Furthermore the major fraction of 
the field pressure is supported by the tension in the field itself and so is called force-free. 
This energy can be accessed by allowing the disks to unwind due to the torque of the 
tension in the field, an electric motor, or more likely because of reconnection dissipated 
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by E • J or rjj 2 . 



3.2 Formation of Astrophysical Force-Free Fields and Their Distribution 

In astrophysics the winding of accretion disks forming nearly every compact object plays 
the role of the conducting disks. The original angular momentum of matter, both baryonic 
and dark is formed randomly by the three body interactions in the first non-linear grav- 
itational collapse starting from initial small perturbations (Peebles 1969). The specific 
angular momentum is large, ~ x 10 7 greater than the limiting specific angular momentum 
of say the MBH of every galaxy. Subsequently a mass selection, Mmbh — W~ 3 M ga i, 
and angular momentum transport process, the Rossby vortex mechanism (Li et al. 2001) 
are invoked to explain how such an accretion disk could form (Colgate et al. 2003). A 
similar angular momentum is formed in pre-stellar structure formation from molecular 
clouds, cores, in star formation, but here the ratio of the specific angular momentum of 
the core to the limiting specific angular momentum of the star is less, ~ 300. Regardless, 
once an accretion disk is initiated, half the gravitational energy of forming the collapsed 
object must be released in the transport of angular momentum of the matter accreted to 
the condensed object. By far the largest free energy is released in forming the SMBH, 
10 8 M & c 2 = 10 62 ergs versus 10 11 M & e s t ar s = 10 59 ergs where e s j ars ~ 10 15 ergs/g, the 
specific binding energy of the average star. Hence we concentrate on the the MBH case 
where we see a direct mechanism for converting this free energy into f-f magnetic field 
energy. 

Fig. 2 shows the sequence of all three processes leading to the generation of the inter 
galactic f-f fields. The pre-collapse from the Lyman-a cloud to the flat rotation curve 
galaxy is not shown, but a mass selection is made at a critical thickness by the Rossby 
vortex mechanism leading to a disk of ~ 10 8 M Q (Colgate et al. 2003) in which angular 
momentum is transported by the co-rotating Rossby vortices shown as the disk with 
vortices (Li et al., 2001). Star-disk collisions produce the helicity to produce an a — uj 
dynamo supplying the polloidal, seed field. The differential winding of the Keplerian disk 
flow produces the f-f helix. The f-f helix extends away from the disk a distance determined 
by winding number, n w ~ 10 s years/(27rr s /c) ~ 10 11 turns, or ~ 10 Mpc. Reconnection 
shortens this distance to radio lobe sizes of a Mpc. 

3.3 Filling the filaments with magnetic flux 

The life time of the radio lobes is roughly the same as the formation time of the MBH, 
~ 10 8 years, and so in a Hubble time, 100 times longer, we expect the fields to evolve and 
fill the IGM quasi uniformly. A simulation of the structure of a local region of the universe 
with magnetic fields has been made by Ryu, Kang, & Bierman (1998) in a 32 x 32 Mpc 
box, shown in Fig. 3. Galaxies with typical galaxy spacings are superimposed. In this 
calculation a primordial magnetic field was assumed and its strength chosen such as to 
allow the evolution of structure and the support of this structure as typically observed in 
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Figure 2: A concept of the formation of massive back holes in galaxies. A galaxy forms as a "flat 
rotation curve" disk, M oc R. At some radius in this mass distribution the column density becomes 
great enough, E cr jt — 100 gem -3 such that the Rossby vortex instability is excited, shown as the 
co-rotating vortices in the disk. The resulting transport of angular momentum allows the interior 
mass to collapse forming the black hole. This simple criterion predicts the mass and mass- velocity 
correlation of MBHs. An au> dynamo forms within this disk because of the differential rotation 
and the hclicity injected by star collisions with the disk. The dynamo, with near infinite gain in 
n t ~ 10 12 revolutions in ~ 10 8 years, provides the (small) polloidal flux out of which the helix is 
formed by the differential winding of the foot prints. This helix is force-free because the field lines 
are primarily axial and azimuthal and so matter, tied to the field lines, falls back to the disk in 
the strong gravity. The helix extends a distance of D ~ n t 2irR g ~ 60 Mpc. Reconnection leads to 
radio lobes at smaller distances. 

filaments, (walls) and voids and with a modulation in density of ~ 20 : 1. The variations in 
the orientation of the fields are due to the instabilities inherent to gravitational structure 
formation. However, what is more important is the bounding pressure of the filaments 
which is derived in the calculation from the infall of matter from the voids. This infall is 
the fundamental instability of structure formation due to gravity and dark matter. The 
dark matter follows its own gravitational perturbations and continually passes back and 
forth through itself. The baryons first follow this gravitational perturbation but then 
can not pass through themselves, and therefore make shocks and pressure. This pressure 
is usually assumed to support the filaments from collapse, but cooling should reduce 
this pressure leading to further collapse of the filaments. Instead, in these calculations 
the magnetic pressure supports the filaments both against further collapse, but more 
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importantly against the pressure of gas continuing to infall from the voids. It is this infall 
pressure that confines the magnetic field. They found that the equilibrium magnetic field 
pressure that retains the necessary density modulation ratio of filament density to void 
density had to be as large as ~ lfj, G. This is the field shown in Fig 3. This is very closely 
the magnetic field strength expected if 10% of the MBH formation energy, 10 61 ergs is 
distributed in a GS box, Vol ~ 10 74 cm 3 . This value of the IGM field is not inconsistent 
with the background Faraday rotation measures of distant polarized sources, provided one 
averages field reversals on a scale much smaller, ~ 1/10 to 1/100 of GS, uses an electron 
density of n e ~ 10~ 5 cm" 3 and includes the volume ratio of filaments to voids. We expect 
the synchrotron glow from these distributed, weak fields will outline the structures of 
the universe when such low frequency arrays as LOFAR are operational. For now, only 
the radio lobes and their minimum energy indicate consistency with this picture of the 
magnetic free energy of the universe. With this view of the dominant free energy in the 
universe in the form of f-f fields, we proceed to an abbreviated theory of reconnection 
acceleration of CRs. 




Figure 3: A cosmological structure simulation 32 x 32 Mpc 3 , by Ryu, Kang, and Bicrman 1998, 
in which magnetic field was tied to the matter and the structure evolved with magnetic pressure 
as well as gravity. The final, near-equilibrium state, was one where the magnetic field of ~ lfi 
G was confined to the filaments by the in-fall pressure of matter from the voids. In our model 
of a magnetized universe the magnetic flux, instead of being assumed and adjusted, is instead 
produced in excess by the winding of the accretion disks and reaches the steady state defined by 
the in-fall pressure and reconnection. Hence we have superimposed on a local region of a filament, 
an approximate distribution of symbolic galaxies. CRs will diffuse between galaxies following lines 
of force until escaping to the voids in ~ 10 8 years. This avoids observation of most of the CRs 
that might have been attenuated by the GZK mechanism. 
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4 Reconnection Acceleration of CRs 



4.1 Dissipation of Force-Free Fields in the Laboratory and Astrophysics 

By way of comparison to shock acceleration, there is a vast literature on the experimental 
observation and its interpretation of dissipation, reconnection of, and particle acceleration 
in f-f fields. The organizing physical principle is the maximization of the dissipation rate 
of the magnetic free energy. The free energy of a f-f magnetic field can be accessed or 
converted to another form, heat or kinetic energy, only by E • J. A resistive origin of 
E = rjj leads to heat, but the accleration of relativistic run-away current carriers of say 
energy 'yi i e m i,eC 2 , leads to a kinetic energy density nji e mi e'fi eC* where the number density 
of current carriers is nj^ fi = J/ ec. In order for the current of run-away carriers to exceed 
the usual or expected much larger slow drift of typical current carrying background plasma, 
and thereby effectively cause E^ ~ 0, some un-identified instability must locally immobilize 
these background plasma charged particles so that due to current carrier starvation a large 
E\\ is maintained, E acce i oc dJ/dt. It is reasonable that the greater rigidity of the run-away 
or accelerated current carriers should circumvent the impedance of such an immobilization 
instability, but it is puzzling why, as we observe in cosmic rays, that the bulk of the particles 
energy centers around c 2 or 7 ~ 1. 

In tokamaks, reconnection sometimes leads to "current interuptions" where all the free 
energy of the magnetic field is transformed into a run-away beam of relativistic ions or 
electrons that sometimes "melts a hole" in the metal vacuum liner. Most often reconnec- 
tion leads to an interchange of flux surfaces and a consequential loss of confinement of 
the hot plasma to the cold walls, a multi billion dollar question for fusion. In stabilized 
pinches or reverse field pinches, more analogous to the above described helical fields and 
to astrophysical fields where the winding number is large, and so where Bj, ~ B z , the flux 
surface topology becomes tangled and random. An experimental demonstration of these 
tangled fields is given in the appendix of Colgate (1978). We suspect, but cannot prove, 
that the maximum dissipation rate occurs when the current carriers are accelerated to a 
velocity, ~ c, or 7 ~ 1. The velocity from further acceleration is then bounded, and some 
optimum rigidity occurs for a given acceleration. High rigidity particles should remain 
in a reconnecting flux filament the longest, thereby gaining the largest energy from the 
£|| of reconnection. The momentum of the current carriers, jrriiC, then inertially carries 
the current following the tangled field lines on average to a larger radius and thereby 
rapidly diffuses the current in radius. This is a very complicated non-linear sequence of 
energy flow and plasma and field deformations, but it contains the necessary phenomena 
of acceleration, diffusion and dissipation. It has not yet been simulated. 

4.2 Particle acceleration by magnetic reconnection 

As mentioned before, radio lobes are clear examples where CRs are being accelerated in 
a magnetized IGM. To understand particle acceleration by reconnection in lobes, it is 
instructive to begin in the resistive MHD limit, even though it is probably not valid for 



11 



radio lobes, given their large size and small resistivity in which case ordinary magnetic 

field diffusion will not be fast enough to account for the magnetic energy conversion. For 
example, in a filament of size of ~ 1 kpc, and resistivity of r] ~ 10 4 cm 2 /s (using an 
electron temperature of 10 6 K), the diffusion time will be L 2 /rj ~ 9 x 10 38 sec, much 
longer than the Hubble time. A very different situation obtains, however, by realizing 
that as the fluids carry the frozen-in fields and move them around, steep field gradients 
could be generated; These result in thin sheet-like current structures, and hence greatly 
reduce the diffusion times. In the Sweet-Parker reconnection picture, in which the current 
layer width is A^ ~ (r^ry) 1 / 2 , where ta ~ L/va, the typical MHD time-scale, and va ~ 

6.6 x 10 s i3 3 .io- 6 / n -6 2 cm/s. Then the rate of energy dissipation is related to the rate of 
convection of magnetic flux into and out of the reconnection region. This time scale (again 
in the Sweet-Parker model) is r sp ~ (r^r^) 1 / 2 ~6x 10 25 sec, but clearly still much too 
long to be relevant to radio lobes. 

The physical conditions of the lobes are, rather, more consistent with the so-called fast 
collisionless reconnection scenario, which has recently been studied in the context of hot 
fusion laboratory plasmas (e.g., tokamaks) and magnetospheric plasmas (e.g., Earth's 
magnetotail). This is because for radio lobes, the ion skin depth (which is understood to 
be closely related to kinetic effects in reconnection), di = c/u> p i ~ 2.3 x 10 10 n_^ 2 cm, is 
actually larger than the resistive Sweet-Parker layer width 

sf L \ 1/2 /3 x l(r 6 G\ 1/2 / n \V4/r7\i/2 

A -~ 21x10 Us) (— — ) (w*) (w) cm - (1> 

In this limit, reconnection is mediated by the kinetic physics to break the flux frozen-in 
condition. It has been suggested and shown that the reconnection rate is then independent 
of the resistivity (e.g., Shay & Drake 1998). The exact dependence of the reconnection rate 
on various parameters (especially di), however, is under debate (Shay et al. 1999; Wang 
et al. 2001; Fitzpatrick 2003). Under a simplified geometry and in 2D, the latter two 
studies suggested that the length of the current layer undergoing reconnection depends 
on the boundary driving, which is unfortunately very difficult to determine in the radio 
lobe situation. Another recent study by Li et al. (2003) on a fully force-free system using 
particle-in-cell simulations has shown that collisionless reconnection that is facilitated by 
the full kinetic physics can indeed proceed at a very fast rate, with flow speeds close to a 
fraction of the Alfven speed. 

The above scale estimates strongly favor the idea that collisionless reconnection in radio 
lobes will be Alfvenic, and could play an important role in converting the magnetic energy 
to particles at a fast rate, given the high Alfven speeds within radio lobes. This may 
therefore be the main mechanism of in situ particle acceleration, as demanded by the radio 
spectral index distributions. More detailed scaling studies, especially the dependence on 
the system size, need to be done before we can model the role of reconnection in radio 
lobes in more detail. 

For the general IGM, a similar comparison between the resistive layer width and the ion 
skin depth can be made. The higher density and likely lower magnetic field strength will 
tend to bring these two scales closer but the key uncertainty will be whether and how 
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thin current sheets can be produced via, say, ideal MHD processes. Many uncertainties 
remain, such as the global magnetic field configurations in the lobes and the general IGM, 
how efficiently thin current sheets can be made, etc. 



5 The Power-law Spectrum 

The accepted theory of cosmic ray acceleration is shock wave acceleration in the ISM 
driven by supernova (Axford, Leer, & Skadron, 1977; Bell, 1978; Blandford & Ostriker, 
1988). "This acceptance has been largely based upon the good agreement between the 
"universal" power-law spectrum predicted by shock acceleration, i.e., the power-law index 
becomes: 

r = (d lnN)/(d InE) ~ -(2 + e) 

depending only on the Mach number and the observed or inferred particle spectra." (see 
Blandford & Eichler, 1987 for a review, and many papers by P. Biermann for a more 
accurate comparison.) This belief that a nearly correct power-law spectral index alone is 
unique is instead, a less restrictive condition than commonly believed. Any accelerator for 
which a fractional gain in energy, d InE, by a few particles is accompanied by a fractional 
loss, — d InN, in number of the remainder will give a power-law: 

dN/N = -T(dE/E) . 

The fractional loss for a fractional gain in energy is what would be expected for a rigidity 
dependent loss mechanism where the probability of a relativistic particle being scattered 
out of an acceleration region is inversely proportional to its energy or rigidity. For values 
of T >~ —2, i.e. a smaller fractional loss, the integral energy becomes asymptotically 
large, and at some energy will truncate or limit the acceleration mechanism, destroying 
the confinement and hence the accelerating mechanism itself. Hence, it is not likely that at 
any one time we should see many such accelerators occurring naturally in the Galaxy. On 
the other hand accelerators with V << — 2 will produce a steep spectrum that, relative to 
another less steep one, i.e. T closer to -2, will be lost relative to the less steep mechanisms 
above some critical energy. Hence it is likely that the spectrum of any observed mechanism 
should be close to T = — (2 + e). 

6 Experiments 

Laboratory experiments can be performed to simulate both magneto-hydrodynamics as 
well as the tearing mode reconnection and the associated En acceleration of the "run- 
away" particles. The spheromak and reverse-field pinch experiments are a step in this 
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direction. Interruptions in tokamaks are already laboratory proof of this acceleration. 
Experimental proof of an a — Q dynamo is similarly needed. We need to perform more 
laboratory plasma experiments to observe reconnection in the collisionless limit. Without 
laboratory experiments, as for example shock acceleration, we are still uncertain about 
the origin of cosmic rays. 
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